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Abstract: We study the phase transition in the holographic entanglement entropy for 
various confining models. This transition occurs for the entanglement entropy of a strip at 
a critical value of the strip width. Our main interest is to examine the critical width for 
models with several parameters. For these models, the critical width, the glueball mass and 
the string tension all become functions of these two parameters. Comparing the behavior of 
the critical width in the entanglement entropy and these other scales, we find that l c seems 
to follow closely the deconfinement temperature and the glueball mass. The behavior of 
the string tension is similar to l c , despite of being parametrically smaller than the other 
quantities. 
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1 Introduction 

Entanglement entropy (EE) Sa is a measure of entanglement between two subsystems A, B 
of a quantum system (for a holographically centered review see [1] ) . It is the von Neumann 
entropy of the reduced density matrix, obtained from the density matrix constructed by 
tracing out the degrees of freedom of our subsystem B (pa = tr^p): 



This quantity is if two subsystems are independent (ie. if the reduced density matrix still 
corresponds to a pure state). One can intuitively think of Sa as the entropy due to lack of 
information of an observer in A who can not access degrees of freedom in B. 

This definition of entanglement entropy is easily extended to quantum field theories 
as follows [2] : if the spacetime dimension of the theory is d + 1 we first need to pick a d 
dimensional time slice or Cauchy surface. Then we define a subsystem A to be a (spatial) 
subregion of this Cauchy surface and denote its complement B. The d — 1 dimensional 
hypersurface which divides these two regions (the boundary of A, OA) is what we will call 
the entangling surface. With these considerations, we can define the entanglement entropy 
of the system as in (1.1). Unfortunately, in a QFT, the result will always be divergent and 
we will have to introduce a cutoff £. The first divergent term is (for theories with d > 1) 
proportional to the area (A) of the boundary of A : 



Sa 



tv A (pa log pa) 



(1.1) 



S A oc 




(1.2) 
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This relation is also called "the area law". However, the coefficient of this term is not 
universal, i.e. it depends on the way we choose the cutoff. After this term there can be 
sub-leading divergent terms (including a logarithmic term) and a finite constant term. The 
coefficient in front of the logarithmic term is universal and for CFT's with d even, it is a 
combination of the central charges appearing in the trace anomaly. This relation between 
central charge and EE has been studied recently in [3-6] . The constant term has important 
uses in the condensed matter community [7, 8]. For topologically ordered systems, where 
there is no notion of order parameter, this quantity has proven to be the observable that 
characterizes the phase of the system. 

The gauge/gravity duality allows us to perform calculations for strongly coupled d + 1 
dimensional field theories using holographic (d + 2 dimensional) gravitational duals [9-11]. 
Ryu and Takayanagi [12, 13] proposed a way to compute the entanglement entropy using 
holography. The holographic entanglement entropy is obtained by considering the area of 
the minimal surface (extended in the extra direction) whose boundary ends in dA 1 . If we 
denote this surface 7 (so that dj = dA), the entanglement entropy will be given by [1] : 

Sa = tttt min .4(7) (1.3) 

This formula has not been proved in general but has been successfully verified with a wide 
range of consistency checks- see for example [14, 15]. Furthermore, the result was recently 
derived for the special case of a spherical entangling surface [16] . Although this definition of 
holographic entanglement entropy was first used for CFTs, we will use the same definition 
for studying nonconformal theories (introducing scales such as a confinement length or a 
flux). It is worth commenting that there may be more than one extremal surface from 
which one must only consider the one with the minimal area. The absolute minima can 
be different if we change parameters of our geometry, so this defines naturally a phase 
transition in the EE, which occurs when two different surfaces have the same minimal 
area. 

Ref [17] studied the behavior of the EE of a strip in holographic confining theories. 
On the gravity side these theories are realized by a compact cycle that shrinks to zero at 
some finite depth in the bulk. The holographic dictionary tells us that the bulk direction 
represents an energy scale and hence the minimal radii where the geometry ends introduces 
a mass gap in the dual theory [18]. As we alluded to above, (1.3) can have more than one 
competing extrema for confining theories. We will obtain a "connected solution" where the 
bulk surface runs between the two boundaries of the strip (and its area will depend on the 
width of the strip), but we will always also have a "disconnected" solution which has two 
independent pieces which go straight from the boundary to the cut-off radius (and whose 
area doesn't depend on the width of the strip). By comparing the areas of the connected 
and disconnected surfaces we will find that they coincide for a particular width of the strip 
l c . That means that if I < l c the EE will be obtained from the connected solution and 
for I > l c from the disconnected. In [17], they then compute the value of l c for different 

technically speaking, it is a saddle point, the area will decrease if we extend it in the time direction. 
However if we first Wick rotate to a Euclidean version of the theory, it is indeed minimal. 
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systems and obtain that l c = 0(1) Aj R . That is, the critical width seems to be closely 
related to the confinement scale. The phase transition in the EE for confining theories was 
also studied for different entangling geometries in [19] with a similar result. This phase 
transition in the EE has been studied in different situations [20-23]. Setting holography 
asides, this phase transition in the EE has also been found in numerical lattice simulations 
of SU(N) gauge theories [24, 25]. 

Motivated by [17], in this note we would like to study what is the relation between the 
phase transition of the EE (connected-disconnected surfaces) and an underlying confinement- 
deconfinement phase transition with multiple parameters. In the holographic models stud- 
ied in [17], there is a single scale and so it is natural to expect that l c ~ 1/A/^ as well 
as l c ~ 1/T C (where T c denotes the deconfining temperature). Our purpose is to test 
these relations in a more general context where the holographic models and hence the dual 
confining theory contains multiple parameters. In this extension, we consider holographic 
models with one extra parameter which may help us understand better this transition. The 
EE phase transition will be determined by a critical width Z c , that is, the width of the strip 
in the boundary when the areas of the two contributions are the same. What we will do 
then is to study the behavior of this width and compare it with that of physical quantities 
that characterize our phase transition. 

The structure of the paper is as follows: In section 2, we set the preliminaries for 
the further calculations, that is, we review the computation of the EE for the strip and 
establish the procedure that we are going to follow in the next section. We also explain 
which quantities we are going to consider to compare with l c . In section 3 we consider 
explicit examples with more than one scale where we study what is the dependence of l c on 
the extra scale, comparing its behavior with that of the other physical quantities presented 
in the previous section. The models that we study can be constructed as the double Wick 
rotation of black hole solutions; the latter would be a charged AdS§ BH, D4-D0 bound 
state and the backreaction of a relevant operator in the perturbative regime. We close the 
paper by discussing the results obtained. 

2 Preliminaries 

In this section we will review how to compute the minimal area that will give us the EE 
for our different holographic systems, for the entangling geometry of the strip. 

The field theories we are going to study will be in a confined state, obtained from the 
finite temperature theories via a double Wick rotation, resulting in a solitonic solution. 
Our field theories will be d + 1 dimensional with a compactified internal dimension whose 
proper length shrinks to zero in the bulk (xd ~ Xd + 2irR). A general metric with the 
desired characteristics will be (we are considering that the boundary is at r = oo): 

ds 2 = -g tt {-dt 2 + dx 2 d _ x + fdx 2 d ) + g rr dr 2 + ds'| (2.1) 

We have included the possible presence of an internal manifold X, that usually appears in 
lOd supergravity solutions (but only appear as internal symmetries in the field theories). 
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Where the compact field theoretical direction closes if we get deep in the bulk (/(ro) = 0). 
We have also included the possible presence of internal dimensions. 

If we are dealing with a lOd supergravity solution with nonvanishing dilaton, the 
prescription for the EE is easier to understand with the Einstein frame, so we have to 
rescale our induced string frame metric such that: 

y5L = e- 2 V^ (2.2) 

Note that we will be integrating over 8 dimensions because if the gravitational theory is 
d + 2 dimensional, our integral will be d dimensional. 

Given a metric with rotational symmetry in the spatial directions (as in our case for 
the d — 1 spatial directions), the computation of the EE is straightforward for certain 
cases. In our calculations we will consider the region A to be a strip and so the entangling 
surface are two flat planes separated by a distance I. Let's consider a strip of size D in 
the X2, X(i-\ directions and / << D in the xi direction, it is clear that (for the minimal 
surface embracing the strip in the boundary) only xi will depend on the r coordinate. 
From these considerations it is easy to write the integral for the area: 



A = J x J d d xe-^^ = D d - 2 2nR jf J d Xl e^ ^ f^\g tt f^ 1 + ^ = 

A/2 , 

=D d - 2 2TrRVol(X) / dxxhxJl + h\r 2 (2.3) 

J -1/2 V 

In this equation r = J^-. In the last step, we factorized the volume of the internal manifold 

X and just consider the integral of functions of r, remembering that if the dilaton is nonzero 

we have to include it in h\. Our task now is to find the minimal surface, so we will only 

have to minimize the integral. 

Because the integral doesn't depend on X \ we can use ordinary functional methods to 

dC 

argue that the "energy" is conserved (E = —r — C) and then obtain the maximum depth 
in the r direction of our surface, if we denote this value of the coordinate r* (h* = h(r*)), 
we can rearrange our integral in one that we know how to compute given the maximum 
depth: 



l 2 

1 _ r°° dr 

2 ~ L T 



2 



A c = D d-2 27r RVol(X) / dr 1 2 = = / drT c (2.4) 

Jr* I / ft* \ 2 Jr* 
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This was the computation for the connected surface, there is also a disconnected surface 
(which does not depend on I), that can be easily obtained from (2.3) by demanding r^ 1 = 0: 

poo poo 

jd = D d - 2 2irKVol(X) / drh 1 h 2 = / drT d (2.5) 

J ro Jr 

Now we are interested in compute the difference in EE between the connected and 
disconnected pieces. Because the two areas diverge in the same way (as the integral is 
extended to the boundary), we can get the difference in the finite piece just by subtracting 
inside the integral: 

poo pr* 

AA(r*) = {F c - F d ) - F d (2.6) 

J r* J tq 

In this way, we will define l c as the width where the disconnected and connected piece 
have the same areas: 

AA{r c ) = 0^l c = l(r c ) (2.7) 
2.1 Physical scales in confining theories 

The physical scales that we will consider for the confining models will be the deconfinement 
temperature (T c ), the glueball mass (M g b) and the tension of the flux string (T s t r ). We are 
interested in comparing these energy scales to l~ . 

• Deconfinement temperature: 

We will have two gravitational backgrounds, corresponding to the confining and de- 
confining phases - both of them are periodic in (at least) one spatial direction and 
the imaginary time direction. Their Euclidean action will depend on the radius of 
confinement (R), the temperature (T) and the external parameter (Y). Following the 
standard approach of Euclidean quantum gravity, the phase with the largest (nega- 
tive) action dominates. Hence equating the gravitational actions at the same R, T 
will give us the deconfinement temperature as a function of R and Y: 

AI(R,T c ,Y) = 0^T c (Y,R) (2.8) 

• Glueball mass: 

We can compute the glueball mass of our theory by considering a massless scalar 
field propagating in a background, see for example [26-28]. Using the ansatz: <j> = 
ij)(r)(3(r)e~' lMt , the scalar dynamics are determined by the wave equation: 

<9m W=gaTdu<k) =o (2.9) 

This wave equation can be easily put into a Schrodinger form, but if we do the crudest 
WKB approximation to get a feeling of how the mass depends on the parameters of 
our theory, we don't really need to do the algebra. In this way, if we consider that 
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the glueball is very massive so we neglect other terms proportional to ip(r), for large 
n, WKB yields [27]: 



M. 



9b - 



-1 




,rr 



,tt 



dr 



(2.10) 



The integral is definite and we have to integrate r from the boundary to the end 
of the geometry. Although M is not the first excitation, it gives us a scale for the 
glueball mass spectrum and because we only worry about what is the scale that the 
glueball is defining, we will reabsorb the n factors in this energy scale {M g b). Abusing 
of language, we will call M g i, the glueball mass, even though it is just the scale of 
the large excitations of the glueball spectra (for large n, that we reabsorbed in its 
definition) . 

• Tension of the string: 

The tension of the string is defined for the confining solution as the effective funda- 
mental string tension redshifted at ro [18]: 



Note however that this quantity will be parametrically different than the previous 
ones because it usually comes multiplied by the 't Hooft coupling A, and the calcu- 
lation we will make will be in the strong coupling regime. 

3 Different models 
3.1 5D flux soliton 

The model is a magnetically charged extension of the first case considered in [17]. We may 
begin with the charged black hole in AdS$, which is a solution to the Einstein-Maxwell 
equations of motion and can be uplifted to type IIB SUGRA (see [29]). Wick rotating this 
solution we get a "magnetically charged" soliton. This model has two scales, the potential 
$ at infinity and the compactification radius of the spatial direction x% . In order to study 
a confinement/deconfinement phase transition, we will construct two solutions from this 
black hole: a soliton with magnetic flux in a compactified coordinate and a black hole 
with constant magnetic flux in the compact direction. We expect that if there is a phase 
transition the confinement phase will be described by the thermal soliton and the black 
hole will represent the deconfined phase. 

As noted, given the solution of the charged black hole, we will do a double Wick 
rotation to get the soliton (in order to keep the potential real, we will reabsorb a imaginary 
unit in the potential) and we will assume periodic imaginary time. The other solution 
will simply be the black hole with constant potential. In this way, for the two cases, the 
coordinates t,x$ should be identified as t — > t + iT~ x ,X3 — > X3 + 2irR. The solutions we 
will be considering are then : 



1 




(2.11) 



2ira' 
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Soliton: ds c = [—dt + dxi + dx 2 + f c dx 3 \ H 2 —t— 

A c = V3*(i- r f 2 )d X3 /c =(i_!f)(i + !f + ^) 

Black Hole: ds 2 d = [-fddt 2 + dx\ + dx\ + + 



A d = V3$dx 3 f d =(l- r -i) (3.1) 



The subscripts c, d which denote confined and deconfined phase respectively and cor- 
respond to the soliton and black hole. Note that, because the X3 direction is compact, the 
gauge potential in the deconfined phase is not pure gauge, we have a nontrivial magnetic 
holonomy. From here it is easy to compute the radius of the compact direction and the 
temperature of the black hole in terms of the parameters of the theory (they have to be 
fixed in order to avoid conical singularities). 

The radius of the compact direction x 3 is: 

1 _ro ( 2+ *L\ 



R L 2 \ r J 

L 2 1 $i? 

= — (1 ± VI - 8x 2 ) witha; = — (3.2) 

Like we said before, there are only two independent parameters (we also introduce a 
temperature implicitly but ro is independent of T). The ones that are physically interesting 
are R and and we have seen that tq(R, <&) has two solutions for every pair of parameters. 
However, because we know that at x = 0, ro = t^j, the only physical solution is , so 
from now on we will omit the ± label. Also, because we want ro to be real, the variable x 
will have to be between and 

The temperature of our system will be determined by identifying periodically imaginary 
time for the black hole: 

T- 1 = — (3.3) 
r d 

Although we have two solutions at finite temperature for the same system, at every 
point T, R, <3? the state will be described by one of the solutions. The physical solution is 
the one which has the least (i.e. more negative) free energy. The difference between the 
free energy of the two solutions goes like: 

AI c _ d oc RT- 1 [rj -4- ($Lr ) 2 ] (3.4) 

If we now plug the expressions for T and R in terms of the parameters it is straight- 
forward to find the deconfinement transition (if T < T c we are in the confined phase): 

h(x) = i(l - 4x 2 - 8x 2 + \A - 8x 2 ) (3.5) 
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Figure 1. h(x) 



1/4 



2irR 



We can see the function /i(x) 1 / 4 plotted in figure 1, because x S [0, 4g] — )• G [1, 
The glueball mass 2.10 and tension 2.11 will be: 



* M ak 1 + ^ — = = /"°° 1 dy (3.6) 

^ 4i? 96V ' A A /(l- 2/ -2)(l+y-2 + x 2 y -4) V ' 

T str = ^(?V = ^(l + x/T^) 2 (3.7) 

27ra' V L / 16i? 2 v y v 7 



where A, the 't Hooft coupling is defined as usual A 



L 2 
2na' 



3.1.1 Stripe EE 

With all the previous things in consideration, the computation of the entanglement entropy 
is straightforward: everything will reduce to a relatively complicated integral of a function 
of r between r* and oo (2.4) , where r* is as usual how deep our surface gets in the r 
direction, and after integrating another function we will obtain the physical width l(r*). 
In order to simplify the final expressions we did some change of variables. The first change 
of variables is to go to a coordinate where the boundary is at and the upper limit of 

the integral is 1. This coordinate is: y = — . We have a family of EE parametrized by $ 

r 

and r* (which is an implicit function of I), to make the computations easier let's assume 

r* = — , the final redefinition that we will make is q = = . In this way, 

M r 1 + VI - 8x 2 

we will be able to factorize the dependence in ro to end up having just a family of integrals 

parametrized by fj, and q, whose values must be between [0, 1] and [0, V2], respectively. 

With all these simplifications, using the usual formula for the EE, the integral that we 

will have to compute to get the entanglement entropy is the following (we set L = 1): 
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S(l,*) = ^A[l,*] 



A[l, *] = r 2 A[l, ?(*)] = r 2 /* ^[y, M (Z), g ] 

JO 



row r y «y 

Note that we explicitly wrote the dependence of the expressions on ro, this factorization 
will make the numerics much easier and comes precisely from the choice of coordinates we 
did. From now on we will drop the dependence on x of q and ro (we know their explicit 
dependence on x,R). With this, we can integrate the functions T (note that the second 
piece of T d in (2.6) has to be integrated between 1 and and ^ between 0, 1 for a set 
of values of \x and obtain the functions AA[l, q], the ones that are of interest. The functions 
to integrate are: 



«r ! _ / 1- (1-^)^4 

' [y ' ql V y V (i - y 6 + (i - g 2 ) y A (-1 + y 2 ) m 4 ) { ' 



y 3 /i 2 




j/V(-l+^)(-l- M 2 -gV) 

(-1 + y2) (-1 + y V) (-1 - yV - <z 2 yV) (-1 - y 2 + y 4 (-1 + (1 - q 2 ) ^)) 

Integrating these functions we can obtain l c {q) by equating AA[l c ,q] = 0. With ev- 
erything we have, now it is relatively easy to get this values. For every x we will have to 
obtain the functions // c (cc) by solving (2.7) numerically. After that, we plug \x c and x in the 
integral for the physical width, so we end up having the function l c (x). We can see these 
quantities plotted in figure 2. Note that although the functions may seem to be diverging, 
they are finite at x = -4=. 

We can see how l c tracks well the other physical three quantities, in particular it 
seems to follow most closely the glueball mass. These results are not conclusive about the 
preferences of l c , because it seems all of the physical quantities of interest to track to order 
one. Of course all of the relevant quantities track each other at the qualitative level. 

3.2 Magnetically charged D4 soliton 

In order to obtain this model we will begin with the soliton M5 brane supergravity solution 
(for example [30]): 

ds 2 n = H~ 1/3 {dx 2 + fdx 2 4 + dx\ x ) + # 2/3 (^ + r 2 dQl) 



3 

L 



j - 1) « H, 1 (3.10) 
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Figure 2. (Colour online) Dependence on x for l c (black) and the scales of our theory (dashed) 
1/T C (black), T~l /2 (red) and M' 1 oc a(x)/ro (blue), the quantities are normalized with their value 
at x — 0, so they all begin at 1. 



where: 



#3 



r I 



(3-11) 



Here the parameter L is specified in terms of the string theory parameters and the near 
horizon limit (L » r) that we used ends up being the usual field theory limit (large iV c 
but finite "t'Hooft coupling"). If we rotate in the 4 and 11 direction with an angle 6 we 
obtain the rotated solution: 



dr 2 



ds 2 u =H 3 lh \dx 2 + fH^dxj) + H 2/d ( — + r 2 dnf) + H 3 1/3 H [dx u + cotantf (H^ 1 - l) dx 4 ] 



H =1 - ~ with q 3 = sin 2 9 (3.12) 
In order to KK compactify the 11th dimension we will decompose it as usual: 

1 2 



ds 2 n = e-^g^dxPdx" + e^ /3 [dx 11 + A„dx^]' 



(3.13) 



So we end up with the result we wanted (we add a constant to the gauge potential to ensure 
A(r o ) = 0): 



ds 2 c = H 3 1/2 \-fH "'"dx'i + H^(-di A + dx{ + dx\ + dx% 

where the dilaton, RR 3-form and RR 1-form are: 

F A = gJ 1 3L 3 dn 4 A 



-2<j> _ n -2rr-3/2rrl/2 



1 



4-i 



H]' 2 Hi' 2 {^ f+ r 2 d^\ 
(3.14) 



/ / 

— dx4 = <&— — dx^ (3.15) 
Ho Hq 
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This solution has a magnetically charged RR 1-form, which is sourced by D6 branes, 
so it is natural to think of this solution as a D4-D6 brane bound state. Physically, the 
parameters we are interested in are the confining radius and the potential at infinity $ (in 
terms of the angle of the rotation $ = tan0). We can easily compute the confining radius 
of this soliton: 

2 2 9T3/2 

Because the physical quantities are $ and R, we can express ro and q in terms of these 
variables : 



4L 3 

4L 3 1 4L 3 $ 2 / 3 



r ° 9R 2 1 + $ 2 



* 9i? 2 (l + $ 2 )(l + $-2)V3 9i? 2 (1 + $2)4/3 ( 3 " 17 ) 

If we now consider the deconfined solution, it'll just be the black hole solution of the 
DA background, with a constant gauge potential and the X4 direction compactified in a 
circle. It'll simply be then: 

ds 2 d = H 3 1/2 [-f(r d )dt 2 + dx\ + dx\ + dx\ + dx 2 4 ] + Hl /2 (jj^ + r 2 dn 2 ^j (3.18) 
Where the dilaton, RR 3-form and RR 1-form are: 

e -2<^ = g- 2 Hl' 2 F 4 = g-^lPdn* A = <S>dx 4 (3.19) 

The temperature of this black hole is: 

3r 1/2 

T =4^V2- ( 3 - 2 °) 
In that way, if we are to compare the action of the two solutions, they have to have 
the same compact circles (ie. time periodic in imaginary time with period T _1 and X4 
with period 2ttR). If we compare these two solutions, from M theory the two solutions 
only differ by a rotation, and because nothing else than the metric depends on the rotation 
(which only acts through the determinant and the norm of F3), the two actions are the 
same (with different parameters ro, r d ). From this, the deconfinement temperature will be 
determined by ro = r d : 

T c' = ^rrr = 2 * R ^ 1 + * 2 ( 3 - 21 ) 

The glueball mass and the tension of the confined string are: 



(r /r) 3 ) 7T V r J 1 ^/ y Z(\ _ y -3) 

M gb ~ 5.41843T C (3.22) 



1 - T l/2„.-l/2i 1 3i? 



str ~27ra' 



H' /Z H' /Z \ =- — 1WT ^rZ~\RT* (3.23) 
3 \ r =r 2ira'2L 3 u c v ; 
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Where in the last line we got rid of all the numerical factors, and we subsituted the 
t' Hooft coupling in 5d A ~ ^7 (it has length dimensions). Note that in this case, for 
$ = 0, 

^str ~ \RP A ^, so although the power of R is not the same that the one for Z c , 

— 1/2 

we are still comparing l c with T str because it is a fundamental quantity in the confining 
theory with length units and its behavior with q could be similar. 

3.2.1 Stripe EE 

Now consider the EE for a strip, if we change to coordinates z = 1/r and we consider that 
our strip is a 8 dimensional surface that is completely characterized by the profile z(x), the 
induced metric is: 



Because there is no explicit dependence on Hq, the problem is the same one that 
without the -DO branes but with different identifications (ro also depends on <£). 

If we are to compare the two competing extremal surfaces, we will find that the critical 
width is: 



r c x ~ 1.1647 J = 4.8786T C (3.25) 




1 1 /2 

In this system both l c and M~ b track T c (note also that l c M g ^ ~ 1), while Tj r goes 
like T 2 . An interesting observation about this case is that the flux can be as large as we 
want, so we can completely separate 1/T C from R. In this way what the EE phase transition 
seems to tell us is that l c goes like M^ 1 or T" 1 . 

3.3 Relevant perturbation 

In this section we will consider our thermal theory to be deformed by a relevant perturba- 
tion. For simplicity we will limit ourselves to first order in the perturbation (i.e. we work 
with large "temperature"), which will allow us to have analytic functions until relatively 
far. We are going to consider as usual a confinement-deconfinement transition. We will 
follow the notation of [31-33] and most of the results here included (before the EE) can be 
found there. 

The Einstein-Hilbert action that we are considering is (we set the radius of the asymp- 
totic AdSs to be L = 2): 



V(<j>) = -3 + 4m 2 (f) 2 (3.26) 



Note that the normalization of the scalar is not the canonical one (which will be just 
40 2 — > ^4> 2 )- m denotes mass of the scalar and is related to the dimension of the operator 
inserted in the boundary in the usual manner [9, 11]: ml = y/A{A - 4). The physical 
dimensions we consider for the operator are A € (2,4) (note that then m? < 0, \m 2 \ < 1). 



- 12 - 



In the regime of small <f> (its smallness is controlled by 4>q), we can solve the equations of 
motion analytically. The metrics 2 (which are basically 3.1 with <3? = and the backreaction 
of a scalar and, again, we periodically identify t, x^'.t — > t + iT _1 ,X3 — > X3 + 2-kR ) and 
scalar profile are: 



x) 2 



<A(x) = <Ao(2x-x 2 ) A / 4 1 F 2 ^,^;l;(l 

ds 2 c = c 2 2 (-dt 2 + dx\ + dx\ + (1 - x) 2 dx|) + G xx dx 2 
ds 2 d = c 2 (-(l - x) 2 dt 2 + dx\ + dx\ + dxf) + G xx dx 2 

— p^V w (3 - 27) 

Here -A(x) and G xx (see [33]) can be easily obtained from the equations of motion (G xx 
is expanded to leading order in </>q): 



- 1 (z-l)djz 
3 7, (22 - z 2 



t>( \ [\ ^y-y 2 ) 2 
B{z) = J dy4> (y) _ 

xx ~ (2x - x 2 ) 2 " 2x~^ ((1 " X) + (1 " X) J + 3(2x-x 2 ) 2 ~3 

For j4(x) we have the freedom to add a constant to B{z) and we choose it so the 
integral is finite at y = 0. 

Like in the previous cases, t and X3 are identified as t — > t + zT _1 , X3 — >■ X3 + 27T.R. In 
order to compare the two theories we will suppose they have the same 4>q, T, R. In this way, 

— 1/2| 

from the usual definition for temperature of a gravitational system (2irT = C2G XX | 
T and R will be [33] : 

Note that in general 5 C 7^ 5d- Because the deformation is the same in both backgrounds, 
the transition will occur at T _1 = 2ttR (or 5 C = 5d)- 

The tension of the string and the glueball mass will be: 



1 , , 1 

'sir 



W V 9x lXl 9u\ x=1 - 2na A=i - 4 (1 + CR ^2 R 2 

-S^agb&o) (3.29) 



m~u = — r dxj- 



9b W In \l (Z ~ 



2 We again follow the conventions of [33], comparing with the previous ones (2.1), (1 — x) 2 — —fgt 

Cl = \/ — gtt, Gxx = Qrr (3^) 
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Because we only calculate up to first order in 0q, for every physical variable, the quan- 
tities that we will be interested in will be: 



X( 



= 1 + cxcbl (3.30) 



X(0) 

For the two quantities before then: 



c M -i = c R + — ;a gb = a + a^cfc (3.31) 

gb Oo 

where in the second case, we expanded a gb in the perturbation. In this way, in order to get 
c M -i, we will only have to do the integral oa for every A. 

gb 

3.3.1 Stripe EE 

What we want to do now is to compute l c , the width at which the disconnected piece of 
the entanglement entropy has the same area than the connected one. First we construct 
the area functionals for the two cases: 



(1 - x^Gli 2 



1 




(1 — X*)C2 

(i - x)4 

(l-x*)cf , 

F d = {l-x)c 2 2 G l J x 2 (3.32) 

In the usual way, the * denotes the tip of the strip, which will be the upper limit of our 
integration and we can obtain the physical width by integrating I. Because we are working 
in the small (f)o limit, we can only trust the results up the first quadratic term in 4>o, which 
will make everything simpler. We will also have to expand then in x* = xq + (j%5x: 

F = 5 c (Fq + (^o-^A + (PqFsxSx) 

£ = S'^io + <fii A + tfoisxSx) (3.33) 

Here we explicitly extracted the dependence in Sr so all the functions Ti, k only depend 
on x,xo-The shift of x* due to the perturbation is linear and the only term that depends 
explicitly on the dimension of the operator is the one with subscript A (the last term comes 
only from expanding in x*). We are interested in the difference of the constant piece of the 
EE for the two solutions and, because the divergent piece will be the same, the integral to 
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compute will be the following (note that for the term linear in 5x there is no disconnected 
piece and it is convergent): 



^(a-*) = r dx(Tf - Ft) - f dxFf 

JO Jx* 



(3.34) 



If we now consider x c , the point where the EE from the two contributions is the same, 
we can do the same expansion x* = x c + <f>\bx c (Ao(x c ) = by definition) : 

0oG4a(O + A &x {x c )5x c ) = (3.35) 

In this way we obtain 5x c (that only depends on A). In order to get l c we only have 
to integrate || until x c . If we expand l c then and read the contribution quadratic in the 
perturbation then: 



(RS R )-H c (<h = 0) " / C (0 O = 0) n 



k Sx = 7^ ~ V \_ ' (3-36) 

'5a; J- X c 

The two integrals in kg x are divergent around x c , but they diverge with the same power, 
so ksx is just the ratio of the coefficients that multiply the divergence. The coefficient q c 
will have contributions both from the integrals and from 5r and will be the quantity of our 
interest. To compute it we only have to do two integrals: /a and Aa for every A. 

The critical width of the unperturbed system x c will be : x c ~ 0.15406. 

In this way, we have all we need to study how does the behavior of l c compare with 
that of the other physical quantities. In figure 3 we can see plotted the four c of interest 
as a function of A and of \m 2 \ (remember that m is a purely gravitational quantity and T c 
doesn't depend on $q). 

In this case we see how the tension doesn't seem to capture as well the behavior of l c 
as the other two quantities. It's worth mentioning that q c = — cr + b (by b we mean the 
slope of (3.36)), so q c considerably smaller than cr (is the same than the tension) requires 
some tuning. 



4 Discussion 



In section 2 we compared l c with three quantities that characterize the confined state: the 
deconfinement temperature, the tension of the string and the energy scale of the highly 
excited glueball states (which we called glueball mass). Overall, we have seen that although 
the behavior of l c is not completely the same for all the cases, there are some observations 
that apply to the three cases. 

First, as we discussed earlier, T s t r is parametrically distinct from other scales in holo- 
graphic models. But, even if T str is not close to l c , 7^*^=0) is similar to yj^=q)- 

Second, M g b and T c are approximately the same (although not identical) and the way 
we obtained M g ^ is of course approximate. Comparing these parameters with the critical 
width, l c follows very closely both scales, being roughly the same. It seems to have a 
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Figure 3. (Colour online) Dependence on Cj on A and |m 2 | for c; c (black) and the scales of our 
theory (dashed) Ci/x c = (black), —c^/r^: (red) and c M -i (blue) 



slight preference for M g b, something which makes sense since correlations in the theory are 
dominated by its excitations. This preference reflects the intuition that correlations are 
minimal beyond l c ~ jj-^. 
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